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Abstract 

Recent news cast doubts on London Interbank Offered Rate (LIBOR) integrity. 
Given its economic importance and the delay with which authorities realize about this 
situation, we aim to find an objective method in order to detect departures in the 
LIBOR rate that from the expected behavior. We analyze several interest rates time 
series and we detect an anomalous behavior in LIBOR, specially during the period of 
the financial crisis of 2008. Our findings could be consistent with data manipulation. 
Keywords '.Financial crisis, LIBOR manipulation, interest rates, information theory, 
permutation entropy, statistical complexity 
JEL classification:E43,E47,C65 


1 Introduction 

London Interbank Offered Rate (LIBOR) was established in 1986 by the British Banking 
Association (BBA). Since its inception it has been a benchmark interest rate for pricing 
several derivative financial instruments traded on exchanges worldwide. Many financial 
contracts link parties obligations to LIBOR evolution. 

The BBA defines LIBOR as “The rate at which an individual Contributor Panel bank 
could borrow funds, were it to do so by asking for and then accepting inter-bank offers in 
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reasonable market size, just prior to 11:00 London time”. Each member of the Contributor 
Panel (selected banks from BBA) submits its rates every London business day through 
electronic means to Thomson Reuters by 11:10 a.nr. London Time. When all this informa¬ 
tion is gathered, rates are ordered and the second and third quartiles are averaged. This 
average is published as the LIBOR rate for each maturity and currency. 

On May 29, 2008 [[15] published an article on the Wall Street Journal casting doubts 
on the transparency of LIBOR setting, implying that published rates were lower than 
those implied by Credit Default Swaps (CDS). Investigations conducted by several market 
authorities such as US Department of Justice, and the European Commission, detected 
data manipulation and imposed severe fines to banks included in such illegal procedure. 

The importance of a good pricing system is based on its usefulness for making decisions. 
As [9] affirmed “we must look at the price system as such a mechanism for communicating 
information if we want to understand its real function”. 

In recent years, the journalist initiated investigation about LIBOR make economists 
to critical examine the evolution of LIBOR rates and compare it with other market rates. 
In this line, [22] documented the detachment of the LIBOR rate form other market rates 
such as the Overnight Interest Swap (OIS), Effective Federal Fund (EFF), Certificate of 
Deposits (CDs), Credit Default Swaps (CDS), and Repo rates. The hypothesize that the 
reasons for this divergent behavior is due to expectations in future interest rates and in 
counterpart risk. [20] study individual quotes in the LIBOR bank panel and corroborate 
the claim by [15] that LIBOR quotes in the US are not strongly related to other bank 
borrowing cost proxies. In their model, the incentive for misreporting borrowing costs is 
profiting from a portfolio position. Consequently, the misreporting could be upward in one 
currency and downward in another currency, depending on the portfolio exposition. The 
proof of such behavior is detected with the formation of a compact cluster of the different 
panel bank quotes around a given point. 

p] track the daily LIBOR rates over the period 1987 to 2008. In particular, this paper 
analyzes the empirical distribution of the Second Digits (SDs) of the Libor interest rate, 
and compares with the uniform and Benford’s distributions. Taking into account the whole 
period the null hypothesis of that the empirical distribution follows either the uniform or 
Benford’s distribution cannot be rejected. However, if only the period after the subprime 
crisis is taken into account, the null hypothesis is rejected. 

Given its extensive use, the economic consequences of a “manipulated” or at least a 
wrong LIBOR could be very important. It gives wrong signals to financial markets and 
general public. It could induce excessive or restricted borrowing (depending on the sense 
of the movement). Additionally, it gives a wrong signal to policy makers, which could not 
observe a credit crunch if interest rates are low 

The above mentioned papers raise questions regarding the integrity of the rate quota¬ 
tions coming from individual banks. 

The aim of this paper is to propose the use of quantifiers derived from Information 
Theory, in order to detect anomalies in data time series. In particular, we use permutation 
entropy and permutation statistical complexity, which are able to capture and discriminate 
different dynamics in time series. These quantifiers, when represented in a Cartesian plane 
conform the Complexity Entropy Causality Plane. It is as a useful visual tool to detect 
changes in data stochastic and chaotic dynamics and thus, possible data manipulation. 
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Although statistical methods itself cannot establish the existence of manipulation, this tool 
could be useful for surveillance authorities such as central banks or exchange commissions, 
as an early warning device to propose in-depth investigation. 

2 Information theory quantifiers 

We analyze time series using quantifiers based on information theory. In particular we 
compute two metrics: the permutation entropy and the permutation statistical complexity 
of each series. Given a discrete probability distribution P = {}, with i = {1,..., M}, 
Shannon entropy is defined in m as: 

M 

S[P] =-^Pilnpi (1) 

i= 1 

This quantifier estimate the average minimum number of bits needed to encode a string 
of symbols based on the frequency of the symbols. It equals zero if the patterns are 
fully deterministic and reaches its maximum value for a uniform distribution. The use of 
informational entropy to study economic phenomena is not new. In fact, IMH3I2 could 
be considered as the seminal papers on this field. More recently, m f!6] apply entropy 
and multiscale entropy analysis to assess the efficiency crude oil price. Alvarez Ramirez et 
al. (2012) also use entropy methods to quantify the dynamics of informational efficiency 
of US stock market during the last 70 years. 

However, analyzing time series only by means of Shannon entropy could be insufficient. 
As stressed in an entropy measure does not quantify the degree of structure or 

patterns present in a process and it is necessary to introduce into the analysis a measure 
the statistical complexity in order to characterize the organizational properties of a system. 

A family of statistical complexity measures, based on the functional form developed by 
mi, * defined in [13] [TO] as: 

Cjs = Qj [P ,P e] n[P] ( 2 ) 

where 7~L[P] = 5'[-P]/S' max is the normalized Shannon entropy, P is the discrete proba¬ 
bility distribution of the time series under analysis, P e is the uniform distribution and 
Qj[p,p e ] i s the so-called disequilibrium. This disequilibrium is defined in terms of the 
Jensen-Shannon divergence, which quantifies the difference between two probability dis¬ 
tributions. m demonstrates the existence of upper and lower bounds for generalized 
statistical complexity measures such as Cjs ■ Additionally, as highlighted in m, the per¬ 
mutation complexity is not a trivial function of the permutation entropy because they are 
based on two probability distributions. A complete discussion about this measures and 
details about their calculation is in [25]. 

The planar representation of these quantifiers is introduced in efficiency analysis in [25] 
and was successfully used to rank efficiency in stock markets Wb commodity markets [26], 
and to link informational efficiency with sovereign bond ratings [211] . 
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2.1 Estimation of the probability density function 

In order to evaluate this quantifiers, a symbolic technique should be selected in order to 
obtain the appropriate probability distribution function. Following [2511271 [261 [Ilj , we use 
the [2] permutation method, because it is the single symbolization technique that considers 
time causality. Given a time series of length N, this technique requires the selection 
of an embedding dimension (pattern length) ( D ) and an embedding delay r. No model 
assumption is needed because Bandt and Pompe method makes partitions of the time series 
and orders values within each partition. This methodology requires only weak stationarity 
assumptions. The probability distribution arises naturally from pattern counting. Since we 
work with daily data, we select r = 1. The selected pattern length should fulfill N S> Dl 
, in order to obtain reliable quantifiers . In order to detect changes in the generating 
process dynamics, we compute the permutation entropy and the permutation statistical 
complexity using a sliding window. We fix a window length of N = 500 for embedding 
dimension D = 4. The step between each window is 5 = 30. We select this step because is 
approximately one month. We obtain a total of 118 estimation windows. 

If a series is purely random, permutation entropy is maximized and permutation sta¬ 
tistical complexity is minimized. Since we work with normalized quantifiers, the maximum 
efficiency point of the CECP is (1,0). 

We argue that, without manipulation the location of each window in the CECP should 
be stable, or at least should not follow no predictable path. 

3 Data 

Following [22], we select several interest rates of the UK area. In particular, we consider 
the LIBOR, OIS and Repo Benchmark and Sterling OverNight Index Average (SONIA), 
all at 3 month maturity. The data span is from 17/05/1999 until 08/09/2014, with a total 
of 3996 datapoints, except for UK OIS rate that starts on 19/01/2004 (2776 datapoints). 
All data were retrieved from DataStream. 

SONIA is the reference rate for overnight unsecured transactions in the Sterling market, 
calculated as the weighted average rate of all unsecured overnight sterling transactions bro¬ 
kered in London by Wholesale Markets Brokers’ Association (WMBA) members between 
00:00 and 15.15 GMT in a minimum deal size of 25 million GBP. 

Overnight Indexed Swaps are financial derivatives used to hedge against moves in 
overnight interest rates 

Repurchase agreements (Repo) are collateralised lending transactions. BBA Repo 
Rates are calculated by a panel of 11 contributors, using an averaging method like LI¬ 
BOR. 

4 Results 

Using the data described above and with the methodology described in Section El we 
compute the permutation entropy and permutation complexity for each of the series using 
a sliding window. The rationale for using a moving sample is to study the evolution of these 
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quantifiers during the period under examination. Our sliding window contains N = 500 
datapoints, the frequency is daily (r = 1) and the step between each window is 5 = 30 
datapoints. According to this we obtained 118 estimation windows. It is worth mentioning 
that we performed our analysis for embedding dimensions D = 4. However, if we use 
different embedding dimensions, results are similar. We believe that, with a pattern length 
D = 4 and 500 datapoints, we are able to capture the dynamics of the series under analysis, 
obtaining consistent estimates of both T~is and Cjs■ With the use of the sliding window, 
incoming perturbances are taken into account in the permutation information quantifiers. 

As described in Section [lj economic theory establishes that, if prices are set in a com¬ 
petitive market and information is quickly embedded into prices, their time series should 
look as a random walk. Using the financial economics argot defined by [6], we say that in 
this situation the market is informationally efficient in its weak form. This situation means 
that, in the context of the CECP, estimations of 'Hs and Cjs should lie on the bottom 
right corner. If this happens, the market is said to be informationally efficient. Otherwise, 
either the benchmark model is wrong or the market is informationally inefficient. 

In Figure[l]we display the Complexity Entropy Causality Plane of UK data. Each point 
reflects the calculation of permutation entropy and permutation complexity for a period 
of the sliding window. We observe that SONIA, Repo and OIS rates occupy the bottom 
right area of the CECP. This area could be associated with random processes that exhibit 
no or low memory. 

On contrary, LIBOR rate spreads throughout much of the CECP. In order to analyze 
the direction of the movement across time, we make a zoom of the CECP focusing on 
the LIBOR rate (see Figure [2]). The numbers on this figure references each estimation 
period of this rate. We can clearly observe that LIBOR rate performs relatively good (i.e. 
following dynamics compatible with a random walk) until the estimation of period 53. 
This period spans from December 2004 until November 2006, just prior to the upheaval 
of the financial markets due to the subprime crisis. If we pay attention to the temporal 
evolution of the Information Theory quantifiers, we clearly see a path towards inefficiency 
areas (i.e. higher permutation complexity and lower permutation entropy). This dynamics 
is followed almost regularly until the estimation of period 98, which spans from February 
2010 until January 2012. In May 2008 the Wall Street Journal published the seminal 
notice of possible LIBOR manipulation m, and subsequently a series of investigations 
were started in order to scrutinize the behavior of individual banks in their LIBOR quotes 
submissions. After the estimation of period 98, the path is reverted. Our estimations reflect 
a slow but constant improvement in the informational efficiency of the LIBOR time series. 
This improvement seems to be consistent with the uncover of the LIBOR scandal and the 
subsequent reform in the LIBOR setting characteristics. Our methodology is not directly 
aimed to detect data manipulation. However, by means of it, we observe a significant 
change in the informational efficiency levels in the LIBOR time series, that overlaps with 
the landmarks of the beginning of the LIBOR manipulation and LIBOR scandal detection. 
In fact, several banks recognized their guilty and paid enormous fines for abusing of their 
market positions as published by [3] and m- Putting this things together, we believe that 
our methodology is able to detect situations in which exogenous forces intervene in the 
formation of prices. 

This is specially important because the inefficiency path has not been followed by other 
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interest rates. If LIBOR behavior is clearly different from other interest rate behavior, what 
is the difference between LIBOR rate and the other selected interest rates? Probably, the 
reason could reside in the very intrinsic characteristics of the LIBOR rate. In fact the 
methodology of LIBOR and the other rates differ substantially. Whereas, OIS and SONIA 
reflect an average of actual transactions, LIBOR rate reflects an self-estimation of the 
borrowing cost of a selected group of banks and do not necessarily reflect any transaction. 

In order to clarify our claim, Figure [3] display the temporal evolution of the permutation 
entropy of all the interest rates. This figure clearly shows that the permutation entropy of 
the rates that reflect actual transactions is relatively high and their movement is bounded 
between 0.8 and 1. On contrary LIBOR rate exhibits a progressive deterioration in its 
permutation entropy, achieving a absolute minimum (Hs ~ 0.39) in period 98. After this 
point there is a steady recovery of higher levels of permutation entropy. 


D=4, t=1, S=30 



Figure 1: Complexity Entropy Causality Plane for D = 4,r = l,d = 30 of UK interest rates 


5 Conclusions 

According to our results, Information Theory quantifiers are suitable instruments for study¬ 
ing financial market stochastic dynamics, and eventually to detect anomalies. In particular, 
we found that unilateral interest rates (LIBOR and Repo) exhibited a dramatic change in 
their underlying generating process, whereas bilateral market based rates (OIS and SO¬ 
NIA) were more or less constant in their stochastic characteristics. We believe that further 
investigation in this topic is desirable. We will continue working in the fine-tuning of our 
model in order to detect a possible structural break in the LIBOR generating process. 
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D=4,t= 1,8=30 
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Figure 2: Detail of the Complexity Entropy Causality Plane for D = 4,r = 1,(5 = 30 of UK 
interest rates, focusing in the movement of LIBOR rate. Numbers refers to the estimation 
windows. 


Entropy evolution of UK interest rates. D=4, t=1, 8=30 



Figure 3: Permutation entropy evolution of the different UK interest rates, computed for 
D = 4,r = 1,(5 = 30 
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